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Parawing vehicles may have urms-ual values of nany of the mass and 
aerodynamic factors affecting dynamic lateral stability and control. 
These unusual characteristics are due in large part to the fact that 
the center of gravity of parawing vehicles is located far below the 
parawing, whereas conventional aircraft usually have the vertical 
center-of-gravity location near the plane of the wings. The present 
thesis is an analytical investigation of the dynamic lateral stability 
and control of a typical parawing vehicle. The analysis was made using 
three-degree-of -freedom, rigid body equations of motion. Stability 
derivatives used in the calculations were obtained from static and 
dynamic force tests of a parawing model with rigid leading-edge and 
keel members. The analysis is treated mainly in terms of the effects 
of vertical center-of-gravity position, since this was found to be the 
most significant factor affecting the lateral stability and control of 



the lypothetical vehicle. 
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V . ITflROnJCTION 

A wide range of applications of the parawing concept is currently- 
being considered. This type of wing can be packaged and deployed in a 
manner similar to recovery parachutes, and in addition provides attract- 
ive glide and landing flare capabilities. Seme proposed uses of the 
concept are cargo dropping, spacecraft recovery, manned utility vehicles 
anii rocket booster recovery. A considerable amount of aerodynamic data 
have been obtained in support of these applications, however, the factors 
affecting the stability and control of such vehicles have not been 
clearly defined. Most investigations of the dynamic stability and 
control of parawing vehicles have been of a qualitative nature -using 
dynamically scaled free-f light models. (See refs. 1 to 4, for exanple.) 
In view of the lack of detailed analysis concerning the stabili-fcy and 
control of parawing configurations, the present investigation was made 
to improve the basic understanding of the major factors influencing the 
stability of such vehicles. The stability and control characteristics 
of parawing vehicles might be expected to differ from those of 
conventional aircraft because certain physical characteristics of the 
configurations differ markedly from the characteristics of conventional 
aircraft. Among these physical characteristics are: (l) a center-of- 

gravi-fcy location far below the wing, (2) a mass distribution in which 
the mass is distributed mainly along a vertical axis, and (j) a wing 
having a shape which tends to produce large values of lateral force due 
to sideslip (Cyp)» The second factor is a result of the first, and the 
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third factor is inrportaiit because of the first. Since the low center- 
of-gravity position was expected to have a predominant influence on the 
dynamic lateral stability and control of a parawing vehicle, the analysis 
is treated mainly in terms of the effect of vertical center-of-gravity 
location on lateral stability. The effects of individual stability 
derivatives were also investigated to determine which ones were primarily 
responsible for the effects of center-of-gravity location on dynamic 
lateral stability. 

The investigation consisted of a theoretical determination of the 
dynamic lateral characteristics of an unpowered parawing-payload 
combination gliding at maximum lift-drag ratio for several vertical 
center-of-gravity locations. Stability derivatives used in the calcu- 
lations were based on the restxlts of static and dynamic force tests of 
a parawing model having an aspect ratio of 2.85 and deployed leading 
edge sweep angle of 50 °* The effects of center-of-gravity location on 
the stability derivatives are considered important in themselves to the 
understanding of the dynamic stability and control of parawing vehicles. 
Consequently, the effects of center-of-gravity location on the deriv- 
atives are analyzed as a distinct part of the investigation. 



9 


VI. SYMBOLS AND NOMENCLATURE 

The calciilated stability and control resiilts are presented with 
respect to the stability axis system shown in figure 1. All basic 
force test data of the parawing model are presented for the body axis 
system shown in figure 2. The aerodynamic coefficients are based on 
the flat-pattern characteristics of the parawing (45° leading-edge 
sweep condition) given in table 1. 


A,B,C,D,E 

b 

"c 

% 

Cl 

Ci 

ACj 

Cm 

Cn 

ACn 

Cy 

AOy 

Cl/2 


I>b 

f 


coefficients of lateral stability quartic 

wing span, ft 

mean aerodynamic chord, ft 

drag coefficient, 

lift coefficient, 

rolling moment coefficient, 1^/q^Sb 
incremental rolling moment coefficient 
pitching moment coefficient, 
yawing moment coefficient, M^/q^Sb 
incremental yawing-moment coefficient 
side-force coefficient, 
incremental side-force coefficient 

number of cycles required for the amplitude of a lateral 
oscillation to decrease by a factor of 2 
differential operator (d/ds-^j) 
frequency of oscillation, cps 
drag, lb 
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g 

k 


Aq 




Kn 


Kn 


1 -« 4>+ 1 >% 

u f XU 

side-f carce, 11 

acceleration due to gravity, ft/sec^ 
reduced frequency parameter, ^ 

radius of gyration in roll alout principal longitudinal 
axis, ft 

radius of gyration in yaw about principal vertical axis, ft 

nondimens ional radius of gyration in roll about principal 

longitudinal axis (kv /b) 

Aq 

nondimens ional radius of gyration in yaw about principal 
vertical axis (k^^/b) 

nondimens ional radius of gyration in roll about longitudinal 
stability axis 


2 2 2 2 
sin q 


nondimens ional radius of gyration in yaw about vertical 

stability axis 

__ _ _ _ \ 

sin T) + cos tj I 



^XZ 

nondimens ional product- ( 
((kz/ - Kx/jsin Ti c( 

l/d 

lift-drag ratio 


maximum lift-drag ratio 

m 

mass, slugs 

% 

rolling mcanent, ft-lb 

My 

pitching moment, ft-lb 

Mz 

yawing moment, ft-lb 


O'lXI 
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p 

p 

r 

R 

S 


t 


V 


x,y,z 

X z 
b 


£ 

’ b 


a 

P 

P 

Pw 

7 


period of oscillation, sec 

rolling velocity, rad/sec 

free-stream dynamic pressure, Ib/sq ft 

yawing velocity, rad/sec 

Routh’s discriminant 

wing area, sq ft 

nondimensional time parameter based on span (Vt/b) 
time, sec 

time required for the amplitude of the lateral oscillation 
to decrease by a factor of 2, sec 
free-stream velocity, ft/sec 
body reference axes 

nondimensional distances along the body reference axes from 
basic moment reference center (see fig. 5) 
nondimensional distances meas-ured frcm 0.5 heel length 
station to center-of-gravity location measured in axis 
system parallel to stability axes (see fig. l) 
any of the mass or aerodynamic tenns in the lateral equations 
of motion 

angle of attach of heel, deg 
angle of sideslip, deg or rad 
rate of change of sideslip, rad/sec 
angle of sideslip of parawing, deg 
flight path angle (see fig. l), deg 
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angle between reference axis and principal axis, positive 
when reference axis is above principal axis at nose, deg 
(see fig. l) 

aaagle of attack of principal longitudinal axis, positive when 
principal axis is abcfve flight path at the nose (fig. l) , 
deg 

root of stability quartic + DX + E = 0 

lateral relative-density factor (m/pSb) 
mass density of air, slugs/cu ft 

ratio of air density at altitude to that at sea level 
angle of bank, deg or rad 

angle of bank of parawing, positive when right wing tip is 
down, deg 

angle of yaw, deg or rad 

ratio of the amplitudes of roll and yaw present in a mode 
of motion 
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dC7 

Cj = — 

aCn 

% 

II 

o 


bCi 

^pb^ 


II 


hCi 

% = ^ 
2V 


^Cy 

^Yr = ^ 



c 

^ 

i'PSi 

II 


In the present investigation the term "in-phase derivative" refers 
to any one of the cscillatory derivatives that is based on the 
coniponents of forces and moments in phase with the angular displacement 
during the oscillatory tests. The term "out-of -phase derivative" refers 
to any one of the stability derivatives that is based on the components 
of forces and moments 90° cjut of phase with angular displacement. The 
oscillatory derivatives of the present investigation were measured in 
the following combinations: 


Cjp sin a 
Cnp sin a 
Cy-p sin a 



Rolling Oscillation Tests 


In-phase 
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CZp + ui* sin a 

Cup + Cnp sin a \ Out-of -phase 
Cy + Cy* sin a 

ip J 


Yawing Oscillation Tests 

p ^ 

/V Jm 

V/ / . \M • Ju. w / • 

tp ir 

Cnp cos a + Cnj. \ In-phase 

Gyp cos a + Oy^, 



cos a 
Cj^ cos a 
Cy- cos a 


\ Out-of -phase 
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VII. AI'IALxblS 
A. Method of Analysis 

The analysis was made for a hypothetical paraglider system 
consisting of a parawing-payload combination having the mass charac- 
teristics of a recently proposed recovery system. The parawing and 
payload were assimed to he rigidly Joined so that there was no relative 
motion between them except for intentional control movements. Longitu- 
dinal trim was assumed to be obtained solely by varying the fore and 
aft position of the center of gravity of the system. The vertical 
c enter- of-gravity location was varied by locating the payload relative 
to the parawing in such a m a nn er as to yield several specified center- 
of-gravity locations perpendicular to the parawing keel member (in terms 
of z/b) while maintaining the proper position parallel to the keel (in 
terms of x/b) required for longitudinal trim. The stability derivatives 
for each center-of-gravity location were then estimated from measured 
derivatives for the wing alone, and calculations were made to determine 
the dynamic lateral stability and control of the vehicle for each 
center-of-gravity location. 

B. Description of Vehicle 

The recovery system was assumed to employ an aspect-ratio 2.83 
parawing having a conical canopy and flat-pattern sweep angle of 45 °. 

In the deployed, or flight condition the parawing fabric was supported 
in a 50° sweep condition by three rigid tubular members which formed 
the keel and leading edges of the parawing. Paxawing and payload weights 
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were too pounds and jtoO pounds, respectively. The payload aerodynamics 
were assiimed to consist only of drag and were considered only in the 
determination of longitudinal trim. The system center-of-gravity 
location was varied in the previously discxissed manner such that three 
different configurations of the original concept were studied. These 
configurations are herein referred to as A, B, and C, (z/h = 0.25, O. 5 O, 
and 0. 75» respectively) . 

The configurations were assumed to use the wing-h ank system for 
lateral control - that is, the wing was hanked about an axis parallel 
to the keel menher to produce the forces and moments required for lateral 
control. This is, in effect, similar to the center-of-gravity shift 
type of control actually used on parawings with suspended i>ayloads. As 
stated in reference when such a system is •used for control the 
Incremental lateral force and mcanent coefficients produced hy wing hank 
may approximately he expressed as 


A 


= Gyp + Cl sin 
ACn = Cnp Pw ^ ^ 

^2 = Czp Pv, + I Cl sin ^ 




(1) 


where sin ^ = sin a sin 

The first terms of the right hand side of the equations arise from 
the fact that a value of sideslip is obtained at the parawing when the 
wing is hanked for keel angles of attack other than zero. The second 
terms are the contributions of the lateral component of the lift vector 
which has been tilted hy hanking the wing. 
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C. Determination of Stability Derivatives 
The lateral stability derivatives of the configurations were 
assumed to be due to the parawing alone. There was virtually no infor- 
mation available on the dynamic stability derivatives, however j so, in 
order to provide reasonable inputs for the dynamic stability analysis, 
an experimental investigation was conducted to determine both the static 
and dynamic stability derivatives of a 0 . 15 ~scale model of the parawing 
assumed for the theoretical analysis. 

A sketch of the model parawing is presented in figure 5- The model 
was constructed of nonporous plastic membrane attached to three equal 
length rigid members such that in the deployed or fli^t condition the 
rigid menbers supported the fabric in a 50° sweep condition. Additional 
information pertaining to the model can be found in reference 5 * 

The static and dynamic force tests were conducted in a low-speed 
tunnel with a 12-foot octogonal test section at the Langley Research 
Center. The tests were made at a dynamic pressure of I .65 pounds per 
square foot which corresponds to a Reynolds number of 0.99 x 10^ based 
on the keel length. The static lateral tests were made for an angle-of- 
sideslip range of ±5°. Forced oscillation tests were made in yaw and 
roll with angular amplitudes of ± 5 ° and for a value of the reduced 
frequency parameter k of O. 25 . A detailed description of the dynamic 
force test equipnnent and the method of obtaining the dynamic parameters 
is presented in reference 6 . The mcanents were all measured about the 
basic moment center shown in figure the static and dynamic 

lateral stability derivatives of the con 5 )lete configurations were 
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obtained "by transferring these measxired wing-alone model data to the 
required center-of-gravity locations hy means of the equations of 
reference 7 which are included in the appendix of this pai>er for 
convenience. 

D. Method and Scope of Calculations 

Calculations were made to determine the dynamic lateral stability 

and control of the three parawing configurations for the condition of 

l/d ^by gliding flight at an altitude of 10,000 feet. 

The lateral stability calculations consisted of the determination 

of the damping and period of the lateral modes of motion. The damping 

of both the oscillatory and aperiodic modes is expressed in terms of the 

damping factor . - — , the reciprocal of the time to damp to one-half 
■'' 1/2 

amplitude. Positive values of this parameter indicate stability 
(positive damp ing) and negative values indicate instability (negative 
damping) . The calculations were made using the equations of motion 
given in appendix B. Stability boundaries in terms of the dihedral 
effect derivative Ci and the directional stability derivative 

“P 

were also determined. The spiral, or static stability boundary is 
described by 

E = 0 (2) 

while the necessary and sufficient conditions for neutral oscillatory 
stability are that the coefficients of the stability quartic satisfy 
Routh’s discriminant set equal to zero 

R = BCD - AD^ - B^E = 0 


( 3 ) 
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and that B and D have the same sign. The significance cf the 
lateral stability boundaries is indicated by a discussion of the modes 
of motion in the plane. 

CalcxLLations were also made to determine which aerodynamic and mass 
parameters affected the lateral stability of the configurations to the 
greatest extent as the center- of-gravity location was varied. These 
calculations were made using the method of reference 8 which is based 
on a Taylor’s series expansion of the roots of the stability quartic 
near the original solution. Although strictly applicable only for 
small incremental changes in derivatives and mass parameters near the 
original solution, the method does afford an insight as to the important 
terms in the equations of motion as the center of gravity and hence the 
stability derivatives and mass distribution are changed. The results 
are presented as slopes indicating the rate of change of the damping 
and frequency of the various modes of the latereil motion with changes 
in any of the parameters in the equations of motion. 

In conjunction with these calculations, the graphical root locus 
technique was used to help visualize the character of the solution of 
the lateral stability quartic as the mass and aerodynamic parameters 
were varied. This method illustrates the path, or locus of the 
solutions of the quartic on the complex plane. Presented in figure 4 
are the features of the ccDq)lex plane as applied to dynamic systems. 
Lines of constant time to halve amplitude (or double amplitude) are 
parallel to the imaginary axis. Radial lines emanating from the origin 
represent constant values of damping ratios and cycles to halve (or 
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double amplitude. Circles with centers at the origin represent lines 
of constant xmdanqped natural frequency. For stability, all roots of 
the quartic must be located in the left half of the ccmplex plane; that 
is, there must be no positive real parts to the solutions. 

The lateral control calculations consisted of calcTilation of the 
lateral motions following a step input of 5° wing-bank angle. The 
equations of motion were solved using numerical integration techniques 
in conjunction with a digital computer. 

All of the foregoing types of calculations were made for configu- 
rations A, B, and C. In an attempt to provide a better understanding 
of the magnitude of the factors causing the effects of vertical center- 
of-gravity position shown by conparison of the restilts for configurations 
A, B, and C, additional calculations were made for configuration A in 
which the values of mass and aerodynamic parameters were changed one- 
at-a-time to the values for configuration B to determine the effects of 
these individual parameter changes on damping and period. 
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VIII. RESULTS MD DISCUSSION 


A. Presentation of Results 

The results of the parawing model force tests are presented in 
figures 5 "to 8. Calculated lateral stability derivatives and mass 


iicr oxix cc 


full-scale configurations are presented in 


figures 9 "to 12 and are tabulated in table 2. The results of the 
dynamic lateral stability calculations are presented in figures 13 
thrcugh l8 with the sensitivity of the dynamic lateral stability of the 
configurations to changes in the stability derivatives and mass param- 
eters being listed in tables 5 and 4. Root locus plots showing the 
changes in the character of the solutions of the lateral stability 
quartic as the stability and mass parameters are changed are presented 
in figures 19 to 28. Time histories of the motions following wing-bank 
input are presented in figures 29 and 30* The two primary results of 
the investigation - the changes in the stability derivatives and mass 
parameters with center- of- gravity location and the effects on dynamic 
lateral stability and control of these parameter changes - are discussed 
separately. 

It should be noted that any theoretical analysis of the dynamics 
of a parawing configuration can be limited by the concept itself. The 
flexibility of the parawing, the means of connection between parawing 
and payload, and the number of degrees of freedom of the system an 
coniplicate and in some cases may render useless any a naly tic attempt 
to predict the dynamics of the system. Consequently, the calculated 
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resulxs contaliied herein do not apply directly to paraving configurations 
or flight conditions other than those investigated. It is believed, 
however, that the major results and trends shown by these calculations 
are indicative of the dynamic characteristics to be expected of vehicles 
which have the following characteristics: very low center- of- gravity 
location, mass distributed mainly along the vertical axis, and a wing 
which has a large value of the lateral-force derivative Cy . These 
characteristics are typical of a wide variety of paragliders, or gliding 
parachutes . 


B. Stability Derivatives and Mass Parameters 
Mass distribution .- The calculated mass and geometric charac- 
teristics of the three configurations are given in table 2. It should 
be noted that for all three configurations the nondlmensional radius of 
gyration in roll (K^j) was larger than the nondlmensional radius of 
gyration in yaw (K^) • This situation is quite different from that of a 
conventional airplane where the radius of gyration in roll is less than 
that in yaw < K^) and leads to near-vertical inclinations of the 
principal axis of least inertia. This difference would be expected to 
lead to substantial differences in the character of the lateral motions 
between parawing vehicles and conventional aircraft. For e xamp le, the 
large values of inclination of the axis of least inertia will lead to 
increased yawing conq)onents in the lateral modes of motion since a 
vehicle will inherently tend to move about an axis of least inertia. 

It whould be noted that in the present investigation, the principal 
longitudinal axis (the axis to which is referenced) was considered. 
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in the conventional manner, to be the axis located in a generally 
forward direction. The only unusual situation was that for the parawing 
vehicles this was the axis of greatest inertia. 

Longitudinal characteristics .- The parawing system was assumed 
(on the basis of the model force test data of figure 5 with an added 
drag increment to account for the payload) to have a maximum L/d value 
of 4.27 at 25 ° keel angle of attack for all center- of-gravity locations. 
The longitudinal stability of the configurations was not investigated 
in this study, but it might be noted that a paraglider configuration 
beccmes more stable longitudinally for increasing values of x/b and 
z/b. 


Static lateral stability derivatives .- The static lateral stability 

characteristics of the parawing model are presented in figure 6 in the 

form of the static stability derivatives Cy^, ^n^^ plotted 

against angle of attack. The iralues of the derivatives were obtained 

from the differences between the values of the coefficients at sideslip 

angles of +5° and -5°. The variation of Cnp (referenced to body axes) 

shows a negative, or destabilizing trend at about 53° angle of attack. 

It should not be inferred frcm this result that a vehicle with this 

parawing would experience a directional divergence at these higher 

angles of attack. It is a well established fact that when an aircraft 

has large negative values of Cj in conjunction with the negative 

P 

values of Cn^ a- directional divergence does not necessarily occur. 

This point is discussed in detail in reference 9« 
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The values of the static lateral stability derivatives (referenced 

to stability sixes) for the three configurations based on the model data 

and the transfer equations of the appendix axe given in figure 9- 

Lowering the center- of- gravity location from configuration A to that of 

configuration C increased the values of both Cno and -C?,,. The 

“P 

increase in the directional stability derivative Cn^ is caused by the 
fact that as vertical center-of-gravity location Cz/b) is increased, 
the horizontal location (5c/b) must also be changed to maintain longi- 
tudinal trim. 


Lc lateral stability derivatives.- The variation of the in- 


phase oscillatory derivatives with angle of attack for the parawing 
model is presented in figure 7* Static data referenced to body axes 
are also presented and are shown to be generally in good agreement 
with the oscillatory data. 

The variation of the out-of -phase derivatives with angle of attack 
is presented in figure 8. The damping- in-roll parameter (0^^+ sin a) 
and the dan 5 )ing- in-yaw parameter (Cn^. - Cn^ cos a) were essentially 
constant over the angle-of-attack range up to the stall, which probably 
indicates that the p portion of the derivatives were negligible in 
the angle-of-attack range below about This apparent lack of p 

derivatives is in contrast with results of forced oscillation tests of 
swept rigid wings (see ref. 10, for example) in which separated flow 
edffects lave led to large values of the p derivatives at high angles 
of attack. The contrasting characteristics of the parawing may have 
been caused by the flexibility of the parawing itself or by leading- 


25 


edge characteristics. In any event, it was assumed that the rate-of- 

change of sideslip (p) derivatives were negligible and the measured 

ccanbinations were treated as pure rolling and yawing derivatives in 

the dynamic stability calculations. For example, the ccnibination 

Cl + Cl’ sin a was assumed to be for purposes of determining 

P p P 

the lateral stability derivatives of the full-scale configuration. 

The calculated effect of center-of -gravity location on the dynamic 
lateral stability derivatives of the three full-scale parawing configu- 
rations is shown in figures 10 and 11. Lowering the center-of-gravity 
location from that of configuration A to configuration C led to a 

positive increase in the yawing moment due to rolling velocity (C„ ) 

T 

and negative increases in the damping in roll derivative (Cj^) and the 
damping in yaw derivative (Cnj.). The changes in the dynamic stability 
derivatives with vertical center-of-gravity location are cpxite large 
because of the ccsnbination of large distance between the center of 
gravity and the wing and large lateral-force derivative Cy^ which 
enter into the transfer of these derivatives as products. In effect, 
increasing the vertical distance between the center of gravity and 
the parawing keel essentially makes the parawing act as a high vertical 
tail (with short tail length) would on a conventional airplane - that 
is, the side-force derivative CJy^ of the parawing tends to be the 
predominant contributor to the derivative transfer equations. 

The incremental values of lateral force and moment coefficients 
produced by 5° wing-bank angle as calculated from equations (l) are 
presented in figure 12. The incremental rolling moment coefficient 
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is increased “by lowering the center- of =gravity position. The 
yawing moment increment is r^gative, or adverse, for all configurations. 


C. Ijynamic Lateral Stability and Control 

The results of the calculations to determine dynamic lateral 
stability are discussed in terins of the stability of the different 
modes of lateral motion. For conventional aircraft, these modes are 
usually a highly damped aperiodic rolling motion known as the roll 
subsidence mode, a li^tly damped aperiodic motion involving yawing 
and rolling known as the spiral mode, and a lateral oscillation known 
as the Dutch roll. Although the lateral modes of the parawing 
configurations may be different in nature in coEDparison with conven- 
tional nomenclature will be iised to identify the various modes. 

Location of roots on the conplex plane .- The nondimens ional roots 
of the lateral stability quartic for each configuration are plotted 
on the ccmplex plane in figure I3. The two aperiodic (real) roots 
are negative and both are therefore stable. Of particular interest 
is the relatively light damping of the roll mode for these configura- 
tions. The ccmplex pairs representing the lateral oscillation are 
near the neutral stability line, indicating low values of danping. 

The damping ratios for configurations A, B, and C are 0 . 0128 , O.0968, 
and 0.166, respectively. 

Stability boundaries .- Lateral dynamic stability boundaries for 


the configvirations are presented as functions of 


figure 1 ^. All configurations are located in a region of the plane 
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where stability exists, as has been previously pointed out by the 
coanplex plane representation. Increasing the directional stability 
derivative Cnp will lead to spiral instability for both configura- 
tions A and B, however, configuration C is spirally stable for all 

positive values of Cn . Increasing -C?. leads to oscillatory 

p P 

instability in all configurations. An inspection of the root locus 
sketches ctf figure 21 shows that in the case of configuration A, the 
Dutch roll roots cross over into the right half of the ccnplex plane 
and become unstable for increasing dihedral effect. The oscillatory 
Instability which occurs in configurations B and C is the result of 
the roll and spiral roots combining to form a long period oscillation. 

Roll subsidence mode .- The calculated danping factor of the roll 
subsidence mode of the configurations is presented in figure I 5 . The 
data indicate that the system became less damped as the center-of- 
gravity location of the system was moved downward from that of 
configuration A to that of configuration C. The dashed line represents 
the classical single degree of freedom approximation to the roll 

subsidence mode damping j 

VI/2 

reasonably accurate approximation for ccmventional aircraft. Conparison 
of the two resxilts indicates large differences between the complete 
three degree of freedcm results and the simple single degree of 
freedom results. This result indicates that the damping of this mode 
of motion is dependent on factors other than those included in the 
simple approximation. 


C 7 V 


1.57 Pyj 


which generally gives a 
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The calculated variation In the damping factor of the roll 
subsidence mode with changes in the various parameters in the equations 


of motion are presented in table 3* Positive values of 


^1 ^2 

— together 

OX^ 


vitli posi.*txvG i.ncr6s.sss In Xj^ ( 02 * nsgs.’tivs vsluss of 


H/2 

8xi 





with negative increases in Xj^) lead to positive (stabilizing) increases 
in the damping of this particular mode of motion. The results show that 
the yawing moment derivatives and inertia parameters can appreciably 
affect the damping of this mode. The results of the additional calcula- 
tions in which the parameters in the equations of motion of configuration 
A were changed one at a time to those of configuration B in an effort to 
determine the relative magnitude of the changes in the roll mode damping 
due to changes in mass and aerodynamic parameters are presented in 
table 4 together with results predicted by the slope method. The results 
show that the decrease in damping of the roll mode as the center-of- 
gravity location changed from that of configuration A to configuration B 

was caused primarily by the destabilizing contributions of the derivative 

2 

Cjip and the mass parameter . The Increased Importance of the yawing 
derivatives in the determination of the damping of the roll mode can be 
explained by the data of figure l6 which presents the ratio of roll to 
yaw amplitude in the roll mode. Values of this parameter usually are 
in the range of JO to 100 for conventional aircraft but for these 
parawing configurations with their associated stability parameters, all 
values were below 5 a-nd indicate the mode was not a pure rolling motion 
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center-of-gravity location (configuration C), the motion consisted of 
almost an equal rolling and yawing - which is quite different from the 
case for ccaiventional aircraft where the motion is almost pure rolling. 

Spiral mode .- The calcxilated results presented in figure A indicate 
that the spiral mode became slightly more damped as the center of 
gravity was lowered with respect to the wing. An inspection of table 4 
indicates that the stabilizing effect was primarily due to the 
stabilizing influence of changes in those derivatives and 0^^) 

that usually affect the spiral stability of an aircraft. The spiral 
root can usually be approximated by A = -^ and those terms which 
most affect this ratio are seen to be of great importance in the 
damping of the spiral mode for the parawing configuration - as they 
have traditionally been for conventional airplanes. 

Oscillatory mode .- Presented in figures 17 and l8 are summaries of 

the Dutch roll oscillation characteristics of the three configurations. 

The inverse cyclic damping - • ■ and roll to sideslip parameter 

^ 1/2 

for each configuration are presented together with the recommended 
boundary of reference 11 in figure 17. The presentation of the Dutch 
roll characteristics in this manner is not necessarily meant to imply 
that parawing configurations should conform to the handling quality 
parameters of conventional aircraft, but it is one of the present 
standards available for evaluating the effect of the Dutch roll 
characteristics on flying qualities and is consequently used for 
purposes of comparison. The results are also compared in figure l8 to 
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the older handling quality requirements (expressed, in terms of damping 
and period of the oscillation) of reference 12 inasmich as the present 
day aircraft requirements tend to present pessimistic views of the 
oscillatory characteristics of vehicles having extremely light wing 
loadings (the equivalent side velocity Vg of such vehicles is much 
smaller than that of conventional aircraft). The results indicate that 
increasing the vertical distance between parawing and payload (varying 
center- of- gravity location from configuration A to C) led to a longer 
period, better damped lateral oscillation, and smaller values of roU- 
to-sideslip ratio. Only configuration C satisfied the present day 
handling quality requirements with both conf igurations A and B lacking 
satisfactory Dutch roll damping. The res\ilts presented in figure l8 
show that only condFiguration A did not satisfy the older handling 
quality requirements. 

Tables 5 and 4 indicate that the increase in damping as the center 
of gravity is lowered is primarily caused by the stabilizing contribu- 
tion restilting from the positive increase in C^p being larger than the 
destabilizing contribution of The effects of on the damping 

of the lateral oscillation and aperiodic modes have been investigated 
in the past (refs. 15 and l4-) where it was shown that positive increases 
in the derivative usually increases the damping of the Dutch roll mode 
while decreasing the damping of the roll subsidence mode. The deriva- 
tive Cnp will almost always be increased positively by increasing 
vertical distance between the parawing and center of gravity. Increases 
in the same geometric variables will usually lead to negative values of 
Kxz as explained previously. 
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It shcjuld "be noted that in table 3 Ibe contributions of , 

Cjp, and Cjp to the Dutch roll dariping reverse signs between configu- 
rations A and B. That is, increases in the derivatives in the usual 

stabilizing sense (negative for Cj , positive for Cn and C; ) lead 

P 3 3 

to less stable Dutch roll oscillations for conf igurations B and C. An 
examination of the root locus sketches of figures 20 and 21 reveals 
that the location of the circular symbols (zeros) change between 
configuration A and B. This changes the nature of the solution of the 
stability quartic; for example, in figure 21 the Dutch roll mode rapidly 
becomes unstable as dihedral effect is increased for conf Iguratlon A 
while the mode becomes more stable for configurations B and C. These 
results are primarily due to the effect of on the distribution of 

damping of the system. Similar results are reported in reference 13 
for large positive values of Cnp- As discussed in reference l4, the 
total dam p ing of the system can be expressed as b/a. The partial 
derivative of the ratio with respect to any parameter Xj_ must 


therefore follow the relation 


V 




0.693 b 6xi 


h/2 ] 


\ 


— / 

6xi / V ^^i / 


+ 2 


Spiral 

mode 


Roll 

mode 


~^l/2 

6xi 




Dutch roll 
mode 


Since ^ — for and must be zero (changes in these deriva- 

ox^ P 3 

tives do not change the total damping but merely redistribute the 
damping among the modes) the stabilizing contribution of the derivatives 
to the roll subsidence damping (as in configuration A) is accompanied 
by a decrease in the Dutch roll damping. 
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Theoretical investigations in the past (see ref. 15, for example) 

1 


have shown that the reversal of the sign of 


^ 1/2 




'^3 


between configu- 

' Dutch roll 
mode 


rations A and B could be predicted by the algebraic sign of the 
expression 

% - 

Positive values of this quantity indicated that increases in the 
effective dihedral parameter would lead to larger values of the 

Dutch roll dangling factor. The reversal of the sense of the contribu- 
tion of to the damping of the lateral oscillation as the center- 

of-gravity was lowered with respect to the parawing was therefore due 
to the positive increase in brought about by center- of -gravity 

position. 

The data presented in tables 5 and 4 indicate that the decrease 
in frequency of the Dutch roll oscillation with increasing values of 
x/b was due primarily to the mass parameters and Yi-^ and the 

stability derivatives Cnp^ and 

The effect of increasing the raduis of gyration in yaw as presented 
in figure 28 should be noted inasmuch as the root locus shows the 
Dutch roll mode becomes unstable for Increases in yawing inertia. This 
fact may be important for parawing configurations having large values 
of yawing moment of Inertia, such as parawing-booster combinations. 

Lateral response to wing-bank control .- The responses of the 


configuration to 5° wing-bank control are presented in figure 2^. The 
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results of the three-degree- of-freedom calculations show that the 
initial roll response was about equal for a1 1 configurations; hut, as 


the motion progressed, the configurations with the lower center-of- 
gravity locations (configurations B and C) showed greater response to 
wi ng-h ank control. The enlarged plot of the time history of the initial 
yawing motions shown in figure 50 indicates that increasing the vertical 
center-of-gravlty location decreased the initial adverse yawing motion, 
which is opposite to what would he expected from a static aerodynamic 
standpoint. This result is caused hy the favorable contribution of the 
product- of- inertia factor to the initial yawing acceleration as expressed 
by the relation: 




^n- 


5Se 

K 2 







(^) 


As pointed out in reference l6 and as can be seen from equation (4) , 
negative values of give a favorable yawing moment contribution 

and can thereby significantly affect the yawing characteristics in an 
aileron roll as the center- of -gravity location is lowered from that of 
configuration A to that of configuration C. The product of inertia 
parameter, as discussed earlier, will almost certainly be negative for 
large vailues of z/b. 
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IX. CONCLUSIONS 

The results of a theoretical Investigation of the dynamic lateral 
stability and control of a parawlrjg conf igtjration having rigid leadirig 
edge and keel members may be summarized as follows : 

1. Increasing the vertical distance between the center of gravity 
and the parawing keel while maintaining the horizontal location required 
for longitudinal trim led to a decrease in the damping of the roU- 
subsidence mode and an increase in the damping of the spiral mode. 

Lowering the center of gravity also caused an increase in the damping 
and period of the lateral oscillation and a reduction in the ratio of 
roll to sideslip, 

2. The two parameters that most affected, the damping of the lateral 
oscillation with vertical changes in center of gravity location were the 
yawlng-moment-due-to-rolling derivative Cnp and the nondimensional 
product of inertia factor 

3. There was little effect of vertical center -of -gravity location 
on the initial roll response to wing bank control but as the motion 
progressed the configurations having the lower centers of gravity 
displayed the greater roll response. The effects of vertical center 

of gravity location on the mass distribution of the system produced 
significant changes in the initial adverse yawing motion following 
the application of wing bank control. Negative values of Kyz caused 
by lowering the center of gravity relative to the parawing tended to 
reduce the initial adverse yaw. 
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A. Stability Derivative Transfer Equations 
The distances x and z correspond to those of figure 3* Zero 
subscripts indicate data measured about the basic reference center 
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The nondlmensional lateral equations of motion referred to a 
stability-axes system (fig. 1 ), axe: 

Roll 

+ ’Scz“b^+) ' * I + 1 

Yaw 

'I' KxzDb^^) = CnpP + | Cn^J>b2^ + ^ Cn^D^ilf 

Sideslip 

2^^b(l>bP + = CypP + I CYpDbi^ + + I tan 7 ijrj 

When 0 q 6 is substituted for 0 , ilf^e for \|r, and PqS 
for 3 in the equations written in determinant form, A must be a 
root of the stability equation 


AA^ + BA^ + CA^ + DA + E = 0 (l) 

where 

A = 8p^3(kxV - KxZ^) 

B = -2m^ Cy^ + Kjj Cj^ - SK^ Cj^ - - KyzCn^j 

^ ■*■ I ■" ^XZ^ip^Yp 

" ~ ^npKxz^p - I *^npCl^ + Kx2CnpCyp - 


- 

C? Cn„Cy 
k P **P 

= I - 


- f Cn^Cl.% - r CZ,C„pCY„ H- i C^CjpGjp 

I "a ’'(‘^*P% ■ %'^V) 


The damping and period of a mode of motion in seconds are given, 

respectively, by the equations tT /o = - ^ - ‘^ 2. ^ and P = ^ ^ 

-*-/ ^ c V d V 

where c and d are the real and imaginary parts of the root of the 
stability equation. 
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TABLE 1. GEOMEB?RIC CHAEACTERISTTCS OF THE MODEL 


Aspect ratio 

Flat 2.85 

Deplqyed 2.57 

CCL 

Flat, sq ft 12.27 

Deployed, sq ft 11. 16 

Mean aerodynamic chord 

Flat and deployed, in 33 *55 

Span 

Flat, in 70.71 

Deployed, in 64.28 

Root chord,, in 50.00 

Sweep angle 

Flat, deg 45.00 

Deployed, deg 50.00 


kk 


Il/TAOC* A1VTr\ A ITTDr\T\'VTVTA 'K^r' "D A A ■hynPrpIT'DG 
I'lMkJjJ -ML1J-/ ,rtJ^VV>l-'J.X«rU'JJ_\-' X riA UTJ.*JXJ X xji U-> 


OF THE COKFIGUEATIOMS 



A 

B 

C 

% 

-0.2556 

-0.2556 

-0.2556 


.072165 

.08275 

.09325 

-.14287 

-.208458 

-.27405 

% 

.0701 

-.06107 

-.1923 

Cy 

^r 

. 0051^ 

.02422 

.04529 

Cz 

ip 

-.08946 

-.17846 

-.35479 


.05054 

.04852 

.07725 

Cnp 

.009 

.0486 

.09898 

Crij. 

-.01757 

-.02452 

-.0352 


.01015 

. 04 ii 8 

.079694 

2 

Kz 

.0163 

.01697 

.018045 

^XZ 

.000681 

-.003266 

- . 00972 

^^b 

5.185 

X 1 

5.185 

Cl 

.757 

.757 

.757 

tan 7 

-.2338 

-.2338 

-.2558 

n 

-18°12* 

7 ° 56 « 

9 ° 12 « 

x/b 

-.0636 

-.1348 

-.2060 

z/b 

.25 

.50 

.75 

- 


(All values are presented with respect to stability axes) 




















TABLE 4.- VARIATION OF THE LATERAL STABILITY OF CONFIGURATION A 
WITH CHANGES IN THE STABIUTY PARAMETERS* 


! 
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Figure 3-“ Gecametry of the deployed parawing model. Dimensions are 

given in Inches. 
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Figure 4.- Features of the complex plane. 
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of the static lateral stability derivatives of 
; with angle of attack. Data referenced to 
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Figure 8.- Variation of the out-of -phase oscillatoiy derivatives with angle of attack. 

referenced to body axes. 
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ABC 

Configuration 


Figure 9*- Variation of the static lateral stahility derivatives with 
center- of-gravity position. Data referenced to stahility axes. 






Figure 10.- Variation of the rolling lateral stability derivatives with 

center-of -gravity position. 
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ABC 

Configuration 


Figure 11.- Variation of the yawing lateral stability derivatives with 

center-of-gravity position. 






















Figure 12.- Incremental values of lateral force and moment coefficients 

produced by 5° wing-bank angle. 
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Figure I5.- Variation of the damping factor of the lateral aperiodic 
modes vith center-of-gravity position. 
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Figure l6.- Varla'tion of the roll-yaw ratio of the roll subsidence mode 

with center- of-gravity position. 
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?, sec 


Figure l8.- Oscillatory characteristics conrpared with the older 
military handling quality requirements. 


ACy. positive 



Figure 19.- Root locus sketches for variations of 




n«gitlve 
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Figure 20.- Root locus sketches for variations of 





/^Cnp positive 
ACp. negative 



Figure 22.- Root locus sketches for variations of 
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negative 



Figure 24.- Root locus sketches for variations of 






AC( negatlvs 



Figure 25.- Root locus sketches for variations of Cj 




Figure 26.- Root locus sketches for variations of 






Figure 28.- Root locus sketches for increasing values of 
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Configuration 

A = .00068) 


B =-.00327) 

C =- 00972) 



Time, sec 


Figure 30.- Effect of center-of -gravity location on initial yawing 

motion of lateral response. 
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